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Abstract 

Wc present a measurement of the branching fraction for the semilcptonic B decay 
B'^ D^i^u, where can be either an electron or a muon. Wc find r(i?*^ 
D~^£~u) = (13.79 ± 0.76 ± 2.51) ns~^, and the resulting branching fraction }3{B^ 
D'^i~iy) = (2.13 lb 0.12 ± 0.39)%, where the first error is statistical and the second 
systematic. Wc also investigate the B*^ — D'^i~j? form factor and the implications of 
the result for | Vcb \ ■ From a fit to the differential decay distribution wc obtain the rate 
normalization |146|-fi3(l) = (4.11±0.44ib0.52) x 10~^. Using a theoretical calculation 
of F£){1), the Cabibbo-Kobayashi-Maskawa matrix element \Vcb\ = (4.19 it 0.45 it 
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0.53 ± 0.30) X 10"2 is obtained, where the last error comes from the theoretical 
uncertainty of Fd{1). The results are based on a data sample of 10.2 fb~^ recorded 
at the T(46') resonance with the Belle detector at the KEKB e^e~ collider. 

Key words: CKM matrix, Semileptonic, B decay 
PACS: 12.15.Hh, 13.30.Ce, 13.20.Hw 
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1 Introduction 



In the Standard Model of electroweak interactions, the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix [1] are constrained by uni- 
tarity. Therefore, experimental measurements of the precise values of the CKM 
matrix elements are important to understand the phenomenology of weak in- 
teractions. 

In the framework of Heavy Quark Effective Theory (HQET) [2], the semilep- 
tonic decay D^l^v is amongst the cleanest modes that can be used 

to measure the CKM matrix element \Vch\- The differential decay rate for 
— > D^l'V can be expressed as [3] 



dy 487r3 



{ms^+niD+fml+iy'' - if'Fliy), (1) 



where Gj? is the Fermi coupling constant, and m^o and mD+ are the masses of 
the B^ and mesons respectively. The variable y denotes the inner product 
of the B^ and meson four- velocities, which is related to g^, the mass 
squared of the lepton-neutrino system: 

2mBomD+ 



and FE){y) is the form factor. By extrapolating the measured \ Vcb\FD{y) to the 
point of zero recoil of the D'^ meson, \Vcb\ can be determined using a theoretical 
prediction of -Fd(I)- In this paper p represents the four-momentum vector of 
the particle in subscript, while p denotes the three-momentum vector in the 
center of mass (CM) frame and piab is the three-momentum vector in the 
laboratory frame. 

\Vcb\ measurements made using B Di^u decays are less precise than those 
made using B — > D*£~u, due to the suppressed decay rate near the point of 
zero recoil, substantial feed down from B — > D*i~9 and the large combinatoric 
background from fake D mesons. Nevertheless, the B — > Dl'P measurement 
provides a consistency check, and allows a test of heavy-quark symmetry [4] 
through the precise measurement of the form factor. Furthermore, this decay 
has one experimental advantage over B — > D*i~u as there is no slow pion 
involved. 

The B^ — ^ D~^£~u decay is preferred to B~ D^i^u because it has much less 
feed-down background from excited charm meson states. Good particle identi- 
fication and vertexing capabilities allow effective reduction of the combinatoric 
background. 
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In this paper we report measurements of the branching fraction of the semilep- 
tonic decay D^i~u, the CKM matrix element |Vc6| and the form factor 

F£)[y). The lepton i can be either an electron or a muon, and the use of the 
charge- conjugate mode is implied. 

The data sample used in this analysis was collected with the Belle detector [5] 
at KEKB [6], an asymmetric e+e~ collider. The data sample has an integrated 
luminosity of 10.2 fb~^ and contains 10.8 x 10® BB pairs. Another data sample 
with an integrated luminosity of 0.6 fb~^ was taken at an energy 60 MeV below 
the T(4S') resonance, and is used as a control sample to check the e'^e~ qq 
continuum background determination. 



2 Belle detector 

Belle is a large solid-angle spectrometer based on a 1.5 T superconducting 
solenoid magnet. Charged particle tracking is provided by a silicon vertex de- 
tector (SVD) and a central drift chamber (CDC) that surround the interaction 
region. The SVD consists of three approximately cylindrical layers of double- 
sided silicon strip detectors; one side of each detector measures the z coordi- 
nate and the other the r — (f) coordinate. The CDC has 50 cylindrical layers of 
anode wires; the inner three layers have instrumented cathodes for z coordinate 
measurements. Eighteen of the wire layers are inclined at small angles to pro- 
vide small-angle stereo measurements of z coordinates along the particle tra- 
jectories. The charged particle acceptance covers the laboratory polar angle be- 
tween 6 = 17° and 150°, corresponding to about 92% of the full CM frame solid 
angle. The momentum resolution for charged tracks is determined from cosmic 
rays and e+e" fi+fi- events to be (o-|p.^|/|pi|)^ = {O.QQ19\pt\f + (0.0030)^, 
where pt is the transverse momentum in GeV/c. 

Charged hadron identification is provided by dE/dx measurements in the 
CDC, a mosaic of 1188 aerogel Cerenkov counters (ACC), and a barrel-like 
array of 128 time-of- flight scintillation counters (TOF). The dE/dx measure- 
ments have a resolution for hadron tracks of 6.9% and are useful for / K sep- 
aration for \'piab\ < 0.8 GeV/c and \piab\ > 2.5 GeV/c. The TOF system has 
a time resolution for hadrons of o" ^ 100 ps and provides t:/K separation for 
\piab\ < 1.2 GeV/c. The ACC covers the range 1.2 GeV/c < \piab\ < 3.5 GeV/c, 
and the refractive indices of the ACC elements vary with polar angle to match 
the kinematics of the asymmetric energy collisions. High momentum tagged 
kaons and pions from kinematically selected D*^ — > D^n'^, — > K~'k'^ de- 
cays are used to determine a charged kaon identification efficiency of 88% and 
a misidentification probability of 9%. 

Electromagnetic showering particles are detected in an array of 8736 CsI(Tl) 
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crystals located in the magnetic volume covering the same solid angle as the 
charged particle tracking system. The energy resolution for electromagnetic 
showers is {aE/Ef = (0.013)2 + (0.0007/E)2 + (0.008/^^/^)2, ^Yiqiq E is 
in GeV. Neutral pions are detected via their 7r° — > 77 decay. The 7r° mass 
resolution varies slowly with energy, averaging cr = 4.9 MeV/c^. For a ± So- 
mass selection requirement, the overall detection efficiency for neutral pions 
from BB events (including the effects of geometrical acceptance) is ~ 40%. 

Electron identification in Belle is based on a combination of dE/dx measure- 
ments in the CDC, the response of the ACC, the position and shape of the 
electromagnetic shower, and the ratio E/\piab\ of cluster energy registered in 
the calorimeter and particle momentum. The electron identification efficiency, 
determined by embedding Monte Carlo (MC) [7] tracks in multihadron data, 
is greater than 92% for tracks with \piab\ > 1-0 GeV/c. The hadron misidenti- 
fication probabihty, determined from Kg — > tt'^'tt" decays, is below 0.3%. 

The 1.5 T magnetic field is returned via an iron yoke that is instrumented to 
detect muons and mesons. This detection system consists of alternating 
layers of charged-particle detectors and 4.7 cm thick steel plates. The total 
steel thickness of 65.8 cm plus the material of the inner detector corresponds 
to 4.7 nuclear interaction lengths at normal incidence. The system covers polar 
angles between 6 = 20° and 155°; the overall muon identification efficiency, 
determined by a track embedding study similar to that used for the electron 
case, is greater than 87% for tracks with \piab\ > 1 GeV/c. The correspond- 
ing pion misidentification probability determined from inclusive Kg tt+tt" 
decays is less than 2%. 



3 Event selection and analysis procedure 

3.1 Event selection 

For hadronic event selection, we require that each event have at least 5 well- 
reconstructed charged tracks, a total visible energy of at least 0.15 times the 
CM energy and an event vertex that is consistent with the known interaction 
point. Continuum background events are suppressed by requiring the normal- 
ized second Fox- Wolfram moment [8] to be less than 0.4. 

This analysis is based on the neutrino reconstruction method, which exploits 
the hermeticity of the detector and the near zero value of the neutrino mass. 
This method was originally developed by the CLEO collaboration for the 
measurement of S — > iriu and p(a;)£i/ decays [9]. 
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We extract information on the neutrino from the missing momentum (pmiss) 
and missing energy (i?miss) in each event. In the CM frame, the total mo- 
mentum of the system is zero, and the total energy is the sum of the two 
beam energies (£"563111) • The missing energy, momentum and missing mass are 
calculated as follows: 



-S'miss — 2£/beam (3) 
Pmiss = -Spj, (4) 

"^miss ~ -^miss ~ IPmissI ) (5) 



where the sums are over all reconstructed particles i in the event. 

In calculating the missing energy and momentum, we identify each charged 
particle using the lepton and hadron identification devices. We also add the 
neutral showers recorded in the ECL that are not matched to any charged 
track. Moreover, we require that the shower shape be consistent with that of a 
photon and the deposited energy be greater than 30 MeV for the barrel region 
and 50 MeV for the endcap. 

If the only undetected particle in an event is the neutrino, the missing mass 
should be consistent with zero. Due to the presence of a tail in the missing mass 
distribution originating from missing particles, an asymmetric requirement 
—2.0 GeV^/c^ < M^jgg < 3.0 GeV^/c"' is applied. In principle, this selects 
events where there is only one undetected particle in the event and that missing 
particle is the neutrino. Neutrinos are usually produced in conjunction with 
charged leptons (e or /i). Therefore, the presence of two or more leptons in 
an event implies that the missing energy and momentum cannot be used to 
accurately reconstruct the neutrino from the B'^ D^i^u decay. For this 
reason, we select events with only one identified lepton with \piab\ > 0.8 GeV/c. 
The momentum requirement helps to reduce the signal efficiency loss due 
to lepton misidentification, and also reduces the backgrounds from hadrons 
misidentified as leptons and random combinations of leptons and D'^ mesons. 

Since the initial state is charge- neutral, the charges of the reconstructed tracks 
should sum to zero if we detect all the particles in the event and have no 
additional tracks. Hence, we rcqiiirc that the net charge, AQ, in the event 
should be close to zero to reject events with other missing charged particles. 
We allow AQ = ±1, as well as 0, to maintain high signal efficiency. 

The neutrino reconstruction can be biased by particles that go undetected 

by passing down the beam pipe. To reject events with such missing particles, 
we require \cos6^^^^J < 0.95, where 9^^.^^ is the polar angle of the missing 
momentum with respect to the positron beam direction in the laboratory 
frame. 



8 



candidates are reconstructed in the K^h^ti'^ decay mode. Kaon 

candidates are required to be positively identified by the hadron identification 
devices. In addition, we require kaon candidates not to be positively identified 
as either a lepton or a proton. If a charged particle is not positively identified 
as a lepton, kaon, or proton, we treat it as a pion. The three charged tracks are 
then geometrically fit to a decay vertex, and we reject combinations which 
do not form a consistent vertex. We further require \pd+ \ < 2.5 GeV/c in order 
to suppress continuum background. We select K~tt~^7t~^ combinations where 
the invariant mass is within 20 MeV/c^ of the nominal D"*" mass (Figure 1). 

To reduce the feed-down background from B — * D*^X£~u decays, we reject 
candidates that are consistent with being produced in the decay D*^ — 
D^TT^. We select 7r° candidates by requiring that the two photon invariant 
mass be within 16.5 MeV/c^ of the nominal 7r° mass, where the energy of 
each photon is required to be greater than 20 MeV. We calculate the mass 
difference Mk-tt+tt+ttO — Mk~-k+tt+ and if it is within 2 MeV/c^ of the nominal 
value for a D*^ decay, the candidate is rejected. 

A variable cos 6 b-di is defined as the cosine of the angle between p^o and 
Pd+£-{— Pd+ +Pe-)- It satisfies the following kinematic relation: 

2\PBo\\PD+e-\ 



The signal events are distributed mostly within the physically allowed region 
\cos 9 B-De\ < 1, while the background events extend to a much wider range. 
We require that candidates have |cos^b-d^| < 1- 

Since the resolution of the missing momentum is better than that of the miss- 
ing energy, we take {Ep,pp) = (|pmiss|,Pmiss) as the four-momentum of the neu- 
trino. Combining the energy-momentum four- vectors for the reconstructed D'^ 
meson, the signal lepton and the neutrino, and using the constraint of energy- 
momentum conservation, we obtain the fully reconstructed B decay variables, 
the beam constrained mass and the energy difference AE defined as 

AE — (Ed+ + E(- + Ep) — -Bbeam, (7) 

Mbc = ^jEl,^-\pD++Pt- + ap,\^. (8) 



We select events with AE in the range —0.2 GeV < AE < 1 GeV; the asym- 
metric requirement is to reject feed-down background from B D*~^X£~P 
decays. In the calculation of (Equation 8) we correct pp by the factor 

a = 1 + AE/Ep, (9) 
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which is equivalent to imposing /S.E — 0. 



3.2 Background sources 

The background sources fall into five categories: combinatoric, correlated, un- 
correlated, misidentified lepton and continuum. 

Combinatoric background: The dominant background in this analysis is the 
combinatoric background in the D'^ reconstruction. The amount of combi- 
natoric background is estimated using the events in the sideband regions 
1.80 GeV/c2 < Mk-^+^+ < 1.84 GeV/c^ and 1.90 GeV/c^ < MK-n+n+ < 
1.94 GeV/c^, shown in Figure 1. 

Correlated background: If a D'^ and a lepton have the same parent B, but 
do not come from the signal decay S° — > D'^i~9, the event is classified as 
correlated background. Processes such as B^ — > D*~^i^u, B — > D**i~u and 
B D^*^7ri~i' (nonresonant) contribute to this source. This background is 
estimated by MC simulation. We use the measured form factor [10] and de- 
cay rate [11] for S° — > D*'^£~P which constitutes the majority (89%) of this 
background. For B — > D**£~P and B — > D^*^7r£~i', the models in refs. [12] 
and [13] are assumed, respectively. Our assumptions for these decay rates [14] 
are based on existing measurements [11]. 

Uncorrelated background: The uncorrelated background consists of events with 
a real from the decay of one B meson and a real lepton from the opposite 
B. We also estimate this background using a MC simulation. Due to the 
charge correlation with the the lepton in this background is usually from 
a secondary decay, and is somewhat suppressed by the lepton momentum 
requirement. 

Misidentified lepton background: In some cases, a hadron track is misidenti- 
fied as a lepton. The amount of this background is estimated from data. We 
treat each hadron candidate track as if it were a signal lepton and weight 
its contribution according to the misidentification probability measured using 
kinematically identified hadron tracks in data. 

Continuum background: The e"'"e~ — >■ qq continuum background is estimated 
using MC continuum events. We also check that this estimation is statistically 
consistent with the off-resonance data. 
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3.3 Signal yield and efficiency 



Figure 2 shows the Mbc distribution after all the event selection criteria 
described above. In this plot the points with error bars represent the on- 
resonance data. The background components are also shown. 

We apply a final event selection requirement Mbc > 5.24 GeV/c^ to define the 
signal region. The overall signal efficiency is 2.69%. Table 1 lists the number 
of signal events and the estimated backgrounds in the signal region. 



4 Measurement of |V^6|-Fd(1) 



The uncertainty on the reconstructed value of y is dominated by the error 
on the measurement of the neutrino four-momentum. We denote the recon- 
structed value as y. The resolution is improved by making the correction to 
the momentum measurement by the factor a given in Equation 9. Using a 
MC simulation, the resolution of y is found to be accurately modeled by a 
symmetric Gaussian with a — 0.03. 

Figure 3 shows the y distribution for data and the estimated backgrounds. 
After all backgrounds are subtracted, we perform a fit to the y distribution 
to obtain |K6|-p£i(l) ^-nd the form factor. In the fit, the function is expressed 
as [15] 

y2 = V — im 

1=1 ^Ni +Z-i=l 



where Nf'^ is the yield in the i-th y bin, cr^* is the statistical error of Nf'^., 
Ninn is the number of bins and 

r dV 

Nj = NbbB{D+ ^ X-7r+7r+)Ts-o j dy— (11) 



is the number of decays in the j-th bin implied by the fit parameters. Here, 
N^B is the number of BB pairs in the sample, B{D^ — > K^-k^ti^) is the 
D'^ — > K'n'^n'^ branching fraction [11], r^o is the B^ lifetime [11], and dV /dy 
is given by Equation 1. We assume a branching fraction B{T{AS) — > B^B^) — 
1/2. The efficiency matrix [e-ij) accounts for the reconstruction efficiency and 
the smearing of y due to the detector resolution. The statistical uncertainty 
in the efficiency matrix is represented by aij. We use ten y bins over the range 
1.00 <y< 1.59. 
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We use a general parametrization, 

Fn{y) = Fn{l){l - pUy - 1) + c^(y - 1)^ + 0{y - if}, (12) 



which is common to many analyses [15-17]. First we fit to a linear form factor 
neglecting terms oi 0{y — 1)^ and higher. We find |Vc6|Fd(1) = (3.83 ±0.35) x 
10~^ and pI, = 0.69 ± 0.14 with a correlation coefficient p{\Vcb\FD{l) , p^) ~ 
0.96, where the errors are statistical only. The is 7.6 for 8 degrees of freedom. 
Allowing to vary in the region cd > 0, we find consistent results. 

Boyd et al. [18] and Caprini et al. [19] provide different form factor parametriza- 
tions based on QCD dispersion relations to constrain the form factor. In each 
case a relation between and cd is obtained, leaving only one free param- 
eter, in addition to | \4;,|-Fci(l), to be extracted from the fit. They also relate 
higher order terms to these free parameters. Using the form factor given by 
Boyd et a/., we find |K6|^i?(l) = (4.14 ± 0.47) x lO'^, p\^ = 1.16 ± 0.25 
and Cd — 1.06 ± 0.28. Using the parametrization of Caprini et al., we find 
\Vcb\FD{l) = (4.11 ±0.44) X 10-2, p% = 1.12 ±0.22 and cd = 1.03 ±0.23. The 
fit results are shown in Figures 4 and 5. 

We use the Caprini et al. form factor determined from the fit to integrate 
the differential decay rate, dV/dy, over y and obtain the decay rate V{B^ — > 
D~^£~u) = (13.79 ± 0.76) ns~^. This decay rate leads to a branching fraction 
B{BO D+(i-v) = (2.13 ± 0.12)%, where the error is statistical only. The 
decay rate F is not sensitive to the choice of form factor parametrization. 
Table 2 summarizes the different fit results. 



5 Systematic uncertainty 

The systematic uncertainties are given in Table 3. The dominant uncertainty 
is the imperfection of the detector simulation for the neutrino reconstruction, 
which is determined by varying a number of simulation parameters including 
the track finding efficiency, the track momentum resolution, the fraction of in- 
correctly reconstructed tracks, the photon finding efficiency, the photon energy 
resolution, the charged kaon identification efficiency, the charm semileptonic 
decay fraction and the fraction. We take a quadratic sum of all compo- 
nents. 

Since we rely on MC simulation to estimate the correlated backgrounds, we 
vary the relative fractions of D*iu, D**iv and D^*^7rii' within the constraints 
of the measured exclusive and inclusive semileptonic branching fractions [11]. 
To determine the uncertainty from the S° — > D*'^£~P form factor, we vary 
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the form factor within the uncertainties of the measurement [10] while taking 
into account the correlations among the parameters. 

The uncor related, continuum and misidentified lepton backgrounds are very 
small, hence even assuming large uncertainties in these backgrounds makes a 
negligible difference in our result. 

The remaining contributions to the systematic uncertainty are the uncertainty 
in the vertexing efficiency, and the uncertainties in the lepton finding 
efficiency, number of BB pairs, — > K~t:^t:^ branching fraction and the 
B^ lifetime. 

We have tested the stabihty of the result by varying several of the selection 
criteria; in each case no significant change is observed. Furthermore, the results 
obtained with the e~ channel and with the ^~ channel are consistent. 



6 Summary 

Using the missing energy and momentum to extract kinematic information 
about the undetected neutrino in the B^ — > D^i^u decay, we have measured 
the decay rate, |V^b|F/:)(l) and the form factor parameters for a number of dif- 
ferent parametrizations. These results, with statistical errors, are summarized 
in Table 2. We have evaluated systematic errors of 12.5%, 12.6% and 18.2% 
on |Vc6|-Fd(1)! p'd ^'^d r respectively. Using the Caprini et al. parametrization, 
we find the decay rate of B^ D^i^u, 

r(SO ^ D+£-p) = (13.79 ± 0.76 ± 2.51) ns-\ 

corresponding to the branching fraction, 

B{BO D+riy) = (2.13 ± 0.12 ± 0.39)%, 

and we find the rate normalization, 

\Vcb\FD{l) = (4.11 ± 0.44 ± 0.52) x 10"^. 

Using Fd{1) = 0.98 ± 0.07 [19], we find 

\Vcb\ = (4.19 ± 0.45 ± 0.53 ± 0.30) x 10"^ 

where the last error comes from the theoretical uncertainty of Fd(1). These 
results are consistent with the result of B° — > D*+e~P analysis from Belle [20] 
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and with the existing measurements [15-17,21-23] 



From the results of D*^e v analysis at Belle, the ratio of and 

Fd*(1) and the difference between p\ and measured to be 



Fb*(1) 



1.12 ± 0.12 ± 0.12 (Linear form factor) 

1.16 ± 0.14 ± 0.12 (Caprini et al. form factor). 



-2 -2 

Pd - Pd* 



-0.12 ± 0.18 ± 0.13 (Linear form factor) 
-0.23 ± 0.29 ± 0.20 (Caprini et al. form factor). 



where the first error is statistical and the second is systematic after removing 
the correlated error terms between the two analyses. These results are in 
agreement with theoretical predictions [19,24,25]. 
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Tabic 1 

Numbers of signal and background events in the signal region after all event 
selection requirements. 



Total vield 


2518 ± 50 


Combinatoric background 


983 ± 22 


Correlated background 


349 ± 14 


Uncorrelated background 


35 ±4 


Misidentified lepton background 


9±1 


Continuum background 


43 ±7 


Final signal yield 


1099 ± 57 



Table 2 



Summary of the results of the dV/dy fit. 



Model 


\Vcb\FD{l)/lQ-'' 


pI 


CD 


r(ns-i) 


xVdof 


Linear 


3.83 ± 0.35 


0.69 ± 0.14 





13.76 ± 0.76 


7.6/8 


Curvature 


q oq + 0.46 
O.OO _ 0,35 


0.69 1 


00 + ^-^^ 

u-u" - 0.00 


13.76 ± 0.76 


7.6/7 


Boyd et al. 


4.14 ± 0.47 


1.16 ±0.25 


1.06 ±0.28 


13.78 ± 0.76 


8.4/8 


Caprini et al. 


4.11 ±0.44 


1.12 ±0.22 


1.03 ±0.23 


13.79 ± 0.76 


8.2/8 
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Table 3 

Summary of the relative systematic errors. 



Source of uncertainty 


A\Vci,\Fd{1) (%) 


A/51, (%) 


Ar (%) 


u reconstruction simulation 


10.6 


9.7 


15.5 


Correlated background normalization 


2.4 


4.4 


1.9 


D* form factor 


1.5 


2.8 


0.9 


Other background normalization 


0.6 


1.8 


0.4 




4 7 


o.o 


o.o 


Lepton finding efficiency 


1.5 




3.0 


^BB 


0.5 




1.0 


Br{D+ if-7r+7r+) 


3.3 




6.7 




1.0 




2.1 


Total 


12.5 


12.6 


18.2 




Fig. 1. Mj^-^+^+ distribution after all event selection criteria except the Mj^-^+^+ 
requirement. The signal region is indicated by (a) and the sideband regions by (b). 
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Fig. 2. Mbc distribution and estimated backgrounds. The points with error bars 
show the data. The shaded component is the continuum and fake lepton background, 
the verticahy hatched histogram is the uncorrelated background, the cross-hatched 
histogram is the correlated background, and the diagonal hatched histogram is the 
combinatoric background. The open histogram is the signal MC, normalized to the 
measured decay rate. 
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Fig. 3. y distribution after the Mbc > 5.24 GeV/c^ requirement. The histogram 
shadings are the same as those in Figure 2. 
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Fig. 4. y distribution and fit. The points with error bars are data after background 
subtraction. The dashed histogram is the fit result for the linear form factor. The 
solid histogram is the fit result for the dispersion relation (Caprini et al.) form 
factor. 
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Fig. 5. \Vcb\F]:i(y) as a function of y. The points with error bars are the data; the 
curves are fit results for the hnear form factor (dashed) and the Caprini et al. form 
factor (sohd) . The shaded band indicates the statistical uncertainty for the Caprini 
et al. form factor fit. 
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